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Using galaxy pairs as cosmological tracers 
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The Alcock-Paczynski (AP) effect uses the fact that, when analyzed with the correct geometry, 
we should observe structure that is statistically isotropic in the Universe. For structure undergoing 
cosmological expansion with the background, this constrains the product of the Hubble parameter 
and the angular diameter distance. However, the expansion of the Universe is inhomogeneous and 
local curvature depends on density. We argue that this distorts the AP effect on small scales. After 
analyzing the dynamics of galaxy pairs in the Millennium simulation [l], we find an interplay between 
peculiar velocities, galaxy properties and local density that affects how pairs trace cosmological 
expansion. We find that only low mass, isolated galaxy pairs trace the average expansion with 
a minimum "correction" for peculiar velocities. Other pairs require larger, more cosmology and 
redshift dependent peculiar velocity corrections and, in the small-separation limit of being bound 
in a collapsed system, do not carry cosmological information. 



I. INTRODUCTION 

One of the main efforts of modern cosmology is to de- 
termine what is responsible for the observed accelerated 
expansion of the universe |3] . Alcock and Paczynski 
[4] proposed a cosmological test (hereafter denoted AP) 
based on the assumption of statistical isotropy around 
any comoving location. For regions of space-time that ex- 
pand with the background, observed angles and redshifts 
can be translated into proper distances using the angular 
diameter distance d^^^z) and the reciprocal of the Hubble 
parameter H{z). Requiring isotropy in proper distance, 
after translating from angle and redshift measurements, 
leads to measurements of the product dA{z)H{z). 

Because radial information comes from redshifts, AP 
measurements are traditionally limited by peculiar ve- 
locities, also known as comoving velocities [6] . These 
add to expansion-driven redshifts, leading to apparent 
anisotropic clustering if redshifts are assumed to be com- 
pletely cosmological in origin, even if the correct dA{z) 
and H{z) are used to analyze redshifts. These redshift- 
space distortions (hereafter RSD) are degenerate with 
the AP effect, removing signal ^5112], unless assumptions 
are made such as the Universe following a FLRW met- 
ric f7l. In fact, it is simply standard convention that 
makes us split redshift into cosmological and peculiar 
velocity components: considering that pairs of galaxies 
move due to local space-time curvature shows that the 
expansion rate and the RSD component can be strongly 
correlated. In the extreme case of bound systems, for ex- 
ample, the combined pairwise velocity is not dependent 
on background evolution, i.e. the expansion-driven red- 
shift difference across a pair is exactly canceled by the 
RSD signal (see Appendix [A]). 

Marinoni and Buzzi |5] recently proposed a method to 
derive cosmological constraints from pairs of galaxies for 



which peculiar velocities can be modeled. They provided 
a fitting formula for the observed distribution of veloci- 
ties, which can then be used to help break the AP-RSD 
degeneracy. They assume the normalization of the galaxy 
velocity distribution to be redshift and cosmology inde- 
pendent, whereas a more recent work by Jennings et al. 
[H] questioned this statement using N-body simulations. 
In the work presented here we investigate this further, 
considering how well pairs of galaxies, selected using dif- 
ferent properties, trace the cosmological expansion. 



We use the Millennium simulation [T] to test how the 
pairwise velocity of galaxy pairs may contain information 
about the background expansion of the Universe. We ar- 
gue that the local density in which the pairs are found 
may affect the amount of information these pairs carry on 
cosmology, because each patch of the Universe expands 
in a way that depends on the local density. Our analysis 
suggests that selecting isolated pairs, as considered by 
Marinoni and Buzzi, can result in average pairwise ve- 
locities more in line with the Hubble expansion, i.e. they 
need smaller, less cosmology dependent peculiar veloc- 
ity corrections. We also find a better match if low-mass 
tracers are used. 



The layout of this paper is as follows: in Sec. |TT] we 
briefly review the AP effect. In Sec. |III| we describe the 
Millennium simulation and the two semi-analytic models 
used 
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we 



Guo et al. (TU] and Font et al. llj. In Sec 
present and discuss the results we obtained by analyzing 
all galaxy pairs regardless of their local density, while 
in Sec. |V] we consider only isolated pairs. The effects 
of varying galaxy properties are studied in Sec 
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Sec. |VII| we compare the results from the two different 
semi-analytic galaxy formation models used. We then 
conclude in Sec. I Villi 
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II. THE ALCOCK PACZYNSKI EFFECT 

Consider a distribution of particles expanding with the 
Hubble flow, in the redshift interval {z — z + Az/2) 

and subtended by an angle A0. Assuming a FLRW cos- 
mology, the proper size of the object perpendicular to 
our line of sight is given by 



(1) 



where dA{z) is the angular diameter distance to the ob- 
ject. The size of the object parallel to the line of sight is 
given by 



d. = 



{i + z)H{zy 



(2) 



where Az is the difference in the redshift of objects clos- 
est and furthest away from the observer and H{z) is the 
Hubble parameter at the central redshift of the distribu- 
tion. 

Assuming that the collection of particles statistically 
does not have a preferred direction with respect to one 
line of sight, then (di) = (^2), allowing a statistical cos- 
mological measurement |3], from a sufhcient number of 
pairs, of 



H{z)dA{z) 



{l + z)M 



(3) 



Note that Az, z and /S.9 are all directly observable quan- 
tities. The AP effect, as described above assumes that 
Az as measured only depends on the cosmological ex- 
pansion. In fact, the relative velocity of pairs of particles 
depends on the local curvature of space, so this is not 
necessarily a good approximation. 



III. THE MILLENNIUM SIMULATION AND 
SEMI-ANALYTIC GALAXY MODELS 



The merger trees, vital for galaxy formation modeling, 
were then constructed by linking each subhalo found in a 
given "snapshot" to one and only one descendant in the 
subsequent output time-slice. 

We used the data from two semi-analytic models of 
galaxy formation based on the Millennium simulation: 
Guo et al. |10j and Font et al. [TT]. Most of our analy- 
sis uses the semi-analytic model developed by Guo et al. 
[TU] which is based on the growth of and merging of the 
population of subhaloes. Within this catalog, each EOF 
group hosts a central galaxy which sits in the minimum 
of the potential of the main subhalo. Other galaxies as- 
sociated to the same FOF group may sit at the potential 
minima of smaller subhaloes or may no longer correspond 
to a resolved dark matter substructure, the latter being 
know as 'orphans'. The last two collectives of galaxies 
are referred to as satellites, although in Guo et al. [lOj . 
the physical processes affecting satellite galaxies only be- 
gin to differ from those affecting central galaxies when a 
satellite first enters within the virial radius of the larger 
system. This is the radius of the largest sphere with its 
center at the center of the FOF group and a mean over- 
density exceeding 200 times the critical value. 

For our analysis, we have varied several parameters 
from the galaxy catalogs, namely redshift, mass of the 
subhalo that hosts the galaxy, stellar mass and r-band 
rest-frame magnitude. Unless differently stated, the red- 
shift shown in our plots corresponds to 2: = 0.989. 



IV. ALL GALAXY PAIRS 

In this section we study the average pairwise velocity 
of galaxies regardless of local density and galaxy prop- 
erties. We shall compare our findings with predictions 
from linear theory to examine general trends, and test 
the possibility of using randomly selected galaxy pairs to 
trace cosmological expansion. 



In order to quantify how the dynamics of galaxy pairs 
may be affected by factors like redshift, isolation ra- 
dius, mass of the halo etc., we have considered a pop- 
ulation of galaxies from the Millennium simulation [T]. 
This traces the evolution of 2160'^ dark matter parti- 
cles of mass 1.18 X 10^ Mq from redshift 127 to the 
present day inside a periodic box of side 500h~^Mpc. 
The simulation assumes a ACDM cosmology with pa- 
rameters based on a combined analysis of the 2dFGRS 
[I2] and the first-year WMAP data [T3]. The parameters 
are = 0.25, fib = 0.045, I^a = 0.75, 71 = 1, = 0.9 
and Ho = 73 km s'^ Mpc'^ 

Data on dark matter particles were stored at 64 dif- 
ferent times. At each output time, the post-processing 
pipeline produced a friends-of-friends (FOF) catalog by 
linking particles with a separation less than 0.2 of the 
mean value. The SUBFIND algorithm 14 was then ap- 
plied to each FOF group to identify all the substructures. 



A. Method 

For each galaxy pair, we compute the comoving sepa- 
ration d, the pairwise velocity V12 and its square u^2- We 
define the pairwise velocity as: 



V12 



dd 
~dt 



{vi - V2) ■ {xi - X2) 



(4) 



where t is cosmic time, and x and v are galaxy positions 
and velocities. Note that, following the Millennium sim- 
ulation, we work in coordinates that are comoving with 
the Hubble flow, hence, V12 represents the peculiar, non- 
Hubble component of the pairwise velocity. In the plots 
that follow, we shall always show the —H{z) d curve and 
denote it as "static solution" . We shall also highlight the 
zero line, which in these plots represents the Hubble flow, 
and denote it as "comoving solution" . Any data point 
above the comoving solution represents pairs where the 
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two galaxies are receding from each other faster than the 
Hubble expansion, while below they are moving towards 
each other in comoving coordinates. 

We group pairs in bins according to their separation 
d and, for each bin, compute the average (W12) and the 



variance var(ui2) 



(W12) of the pairwise velocity. 



Our definition of expectation value is: 



Wo 



pairs 



(5) 



where iVpairs is the number of pairs in the separation bin. 
In all the plots in this paper, we shall always show (W12) 
as a function of pair separation, with the error bars for 
each bin taken as A/var(tii2)/A^pairs- 

By default, we employ a logarithmic binning in galaxy 
separation. In order to better visualize the data, however, 
in some plots we combine underpopulated bins together 
so that each bin represents at least a minimum number 
of galaxy pairs, iVmin- In this section we use A'^min > 1000 
while, due the poor statistics, we shall employ A^min > 2 
for some of the "isolated" plots in Sections |V] and |VI[ 



B. Results 

In Fig. [1] we show the average pairwise velocity V12 of 
all the galaxies within the Guo et al. [1^ semi-analytic 
model at redshift z = 0, as a function of separation d. 
The velocity curve is represented by the red dots with 
one-sigma error bars, while the blue and black lines are, 
respectively, the static and comoving solutions. 

We also plot the prediction of linear perturbation the- 
ory for V12 as the green dashed line, obtained using the 
prescription from 1151 116) : 



Vl2{d) 



Tr2 



dk k Pm{k)ji{kd), 



(6) 



where we have chosen a bias of 6 — 1, which is in 
agreement with that measured from clustering within the 
galaxy catalogue. 

For separations larger than d > 10 h^^Mpc, we can see 
that the average peculiar velocity is correctly predicted 
by linear theory. As we follow the velocity curve into the 
nonlinear regime, galaxy pairs approach the static line 
until they cross it at d ~ 3h~^Mpc. This crossing marks 
the beginning of the infall regime, where the galaxies 
in the pairs get closer to each other, but with smaller 
velocities as their separation decreases. On the smallest 
scales, the pairs asymptote to the static solution. 

To use galaxy pairs as tracers of the cosmological ex- 
pansion, we need their peculiar velocity to be small with 
respect to the Hubble flow or modellable. A smaller cor- 
rection is required if V12 is closer to zero comoving veloc- 
ity than to the static solution. As we noted above, only 
for d > 10h~^Mpc are the velocities closer to the comov- 
ing solution than the static solution, which is the regime 
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FIG. 1: Average relative velocity for all galaxy pairs at 
redshift z — 0. The blue solid line represents the static 
solution, followed by pairs that have already virialized 
and do not feel the background Hubble flow. The green 
dashed line represents the linear theory model for pair- 
wise velocities according to the prescription found in Ref- 
erences [m [IS] with a bias of 6 = 1. Error bars, which 
are too small to be clearly seen are plotted at the one 
sigma confidence limit, assuming Poisson statistics. 



of linear perturbation theory. On scales d < 3h~-'^Mpc, 
galaxy pairs follow closely the static solution on average. 
Their peculiar velocity component is equal and opposite 
to the Hubble fiow, therefore they do not carry any cos- 
mological information (refer to Appendix [A|. 



V. ISOLATED PAIRS 

In this section we investigate how the dynamics of 
galaxy pairs changes when an isolation criterion is im- 
posed and how this depends on the isolation radius, the 
allowed number of galaxies within this radius, and red- 



shift. We use the same galaxy sample as in Sec. IV 



A. Method 

We initially define a galaxy pair to be isolated within a 
radius riso if each galaxy in the pair has exactly 1 neigh- 
bor within Tiso, and that neighbor is the other galaxy 
in the pair. This is equivalent to drawing two spheres 
of radius rigo centered on the galaxies, and imposing 
the absence of galaxies extraneous to the pair in each 
of the spheres. We shall weaken this requirement, al- 
lowing for the maximum number of neighbors A^neigh in 
each sphere to be larger than 1. Thus we can interpo- 
late between the dynamics of galaxy pairs in the fully 



isolated case (iVneigh = 1) and in the unconstrained case 
(-^neigh — ^ oo). We implement such isolation criterion 
with arbitrary number of neighbors in a two-step pro- 
cess. We first determine the number of neighbors for 
each galaxy in the simulation, and later use this infor- 
mation to select only those pairs where each galaxy has 
less than A'noigh neighbors. 

We fix the isolation radius to riso — 4h^^Mpc, which 
matches the definition of isolation adopted by Marinoni 
and Buzzi Note that in the plots that follow, we 
only look at separations less than the isolation radius to 
ensure that the pair is truly isolated. 



B. Results 

In the upper panel of Fig. [2] we present the relative 
motion of galaxy pairs as a function of their separation 
for galaxies that are isolated within a 4h^^Mpc radius. 
The galaxies are taken from the Guo et al. [10] semi- 
analytic model applied to the Millennium simulation at 
redshift z = Q. The data points are plotted in red with 
one-sigma error bars. The blue line is the static solution 
and the black line the comoving solution - see Sec. |IV A| 
for their definition. 

Our first comment on Fig. [2] regards the error bars, 
which are much larger than those in the non-isolated case 
of Fig. [ij The reason is that the imposition of an isolation 
criterion results in a drastic reduction of the galaxy pairs 
found, which in the case of Fig. [2] are only 694 This 
number is in line with Marinoni and Buzzi [5], who find 
721 pairs for their low- redshift SDSS sample, and 509 for 
then DEEP 2 sample. 

The most striking feature about the dynamics of iso- 
lated pairs is the roughly logarithmic growth of the pe- 
culiar velocity V12 for scales larger than ^ 0.2h^^Mpc. 
The behavior of V12 can be explained when we recognize 
that the dynamics is determined by the combined effect 
of two competing forces: the mutual attraction between 
the two galaxies, dominant for small separations, and the 
disrupting outflow from the void - the void effect, dom- 
inant for separations close to the isolation radius. For 
small separations, the mutual attraction of the members 
of the pair overcomes the void effect and we see an infall 
regime. As we study objects with larger separations, the 
void effect becomes dominant and we see a logarithmi- 
cally growing pairwise velocity vi2. 

In the lower panel of Fig. [2] we plot V12 / Hd, that is the 
ratio of peculiar velocity to Hubble flow. For separations 
0.4 < d < 4h~^Mpc we find an almost comoving regime 
where the peculiar velocities are less than 20% of the 
Hubble flow. In such a regime, the RSD corrections are 
small, and one could hope that they are easier to model 



Equivalent to one isolated pair every 10*^ other pairs for d = 
Ih-lMpc 
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FIG. 2: Upper panel: average pairwise velocity for iso- 
lated galaxy pairs at redshift 2 = 0. The isolation radius 
is taken to be 4h~^Mpc for each member of the pair. 
The blue solid line represents the static solution, showing 
the virialization of pairs. The three shaded regions repre- 
sent three different regimes. The left blue area represents 
the virialization regime: galaxies within these separations 
have virialized and do not experience the background ex- 
pansion. The middle red area shows the infall regime, 
where galaxies start to collapse to form bound systems. 
The right green region shows what we denote as the 'void 
effect': on average, the isolated pair feels a stronger grav- 
itational pull separating the pair rather than making it 
closer. The error bars shown are the one sigma confi- 
dence limit, assuming Poisson statistics. Lower panel: 
ratio of average pairwise velocities to the static solution 
Hd. Note that the y axis of this panel is in logarithmic 
scale. The error bars shown are the propagated one- 
sigma errors from the upper panel. 



so that V12 becomes a proxy of the expansion rate. In 
the following subsections we shall investigate how the co- 
moving regime depends on the isolation radius, the local 
density and redshift. 



1. Varying the isolation radius 

It is interesting to investigate whether the void effect, 
giving the approximately logarithmic growth of Vi2{d) 
for isolated pairs, is still present when we allow the size 
of the isolation radius to vary. We demonstrate this in 
Fig. [3j where we show the peculiar velocities of galaxy 
pairs with riso varying from 0.1h~^Mpc to 4h~^Mpc, at 
redshift z = 0.989. 

For riso ^ 0.6 h~^Mpc the presence of the void severely 
affects the dynamics of galaxy pairs. The logarithmic 
growth of the peculiar velocity is visible, even though it 
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FIG. 3: Average pairwise velocity for isolated galaxy 
pairs at redshift z — 0.989 for different isolation radii 
ranging from O.lh^^Mpc to 4h^^Mpc. 



is just a hint for the rjso = 0.6h~^Mpc data points. The 
cosmological scale do, defined as the separation where the 
peculiar velocity v\2 vanishes, occurs at do — 0.8h~^Mpc 
and appears to be independent of the isolation radius. 

At small separations, we cannot see any noticeable dif- 
ferences between the various rjso datasets. They all follow 
the static solution hne for d < 0.05 h^^Mpc, suggesting 
that isolated galaxy pairs tend to virialize on the smallest 
scales just as non-isolated ones do. 

The number of isolated galaxy pairs at z = 0.989 in- 
creases from 435 to 71,201 when reducing the isolation 
radius from 4 to 2 h^^Mpc, a factor of roughly 160. As we 
have noted above, the riso = 2h~^Mpc pairs still experi- 
ence a regime where peculiar velocities are negligible with 
respect to the Hubble flow. Hence their velocity differ- 
ence still traces the cosmological expansion, although the 
maximum separation for which this is the case is halved 
with respect to the Hso = 4h~^Mpc case. We conclude 
that using pairs isolated within a 2h~-'^Mpc radius as 
cosmological tracers would drastically reduce the statis- 
tical error with respect to the 4h~^Mpc case, while still 
needing minimal corrections for RSD, provided that the 
cosmological dependence of the RSD could be modeled. 
We shall discuss this in more detail later. 
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FIG. 4: Average pairwise velocity for isolated galaxy 
pairs at redshift z — 0.989 for different number of neigh- 
bors within the isolation sphere of 4h~^Mpc. The blue 
solid line represents the static solution, showing the viri- 
alization of pairs. The green dashed line represents the 
average pairwise velocity for non-isolated pairs. The er- 
ror bars shown are the one sigma confidence limit, as- 
suming Poisson statistics. 



non-isolated galaxy pairs, as already shown in Fig. [T] as 
a dashed green curve. As iVneigh increases, the different 
Vi2 curves monotonically fill the gap between the pure 
isolated case and the non-isolated case. A good agree- 
ment between the dynamics of pairs with and without 
isolation criterion is reached once we allow each galaxy 
in the pair to have 5000 other neighbors. 

In the A'ncigh — 10 case, we found 250, 670 pairs, 
roughly a factor 600 more pairs than in the fully iso- 
lated case of iVncigh = 1- Nonetheless, the Vi2 curve for 
-^ncigh = 10 is strikingly similar to the iV„eigh = 1 one. In 
particular, the void effect still seems to trigger the loga- 
rithmic growth of the peculiar velocity, with vi2 crossing 
the zero line at dp ~ 1 h~^Mpc (in the fully isolated case, 
we have do — 0.8h~^Mpc). Hence, we suggest that pairs 
that are not completely isolated trace the cosmological 
expansion almost as well as the fully isolated ones, with 
the added benefit of a much better statistics. 



2. Varying the isolation density criterion 

Here, we investigate how the dynamics of galaxy pairs 
changes if we relax the isolation criterion by increasing 
the allowed number of galaxies A^ncigh within the isolation 
sphere of 4h~^Mpc. 

In the linear plot of Fig. |4j we present the average pe- 
cuhar velocity vi2 at z = 0.989 for 7 values of A'^noigh 
ranging from A'noigh — 1 (equivalent to the pure isolated 
case of Fig. [2]) to iVncigh = 5000. We also plot vi2 for the 



3. Varying redshift 

We illustrate the redshift dependence of the pecu- 
liar velocity for isolated pairs in the top panel of Fig. 
[5] where we plot Vi2{d) for the four redshifts z — 
0,0.5085,0.989,1.504. It is remarkable that for separa- 
tions d > 0.2h~^Mpc the peculiar velocity depends only 
slightly on redshift. The scale do, defined as the separa- 
tion where v\2 vanishes, ranges from 0.6h~^Mpc at z = 
to 0.9h~^Mpc at z = 1.504. This is a small variation if 
we consider that at z = 1.504 the Universe was at one 
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third of its current age. 
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FIG. 5: Upper panel: variation of the average pairwise 
velocity of isolated galaxy pairs with redshift as a func- 
tion of separation. The isolation radius is taken to be 
4h~^Mpc for each member of the pair. Lower panel: 
range in pairwise velocity for each separation bin over 
range in the static solutions at each separation bin. Note 
that the y axis of this panel is in logarithmic scale. The 
error bars shown are the propagated one-sigma errors 
from the upper panel. 

In the bottom panel of Fig. [5]we plot the ratio between 
the range in Vi2 and the range in Hd at a given separa- 
tion. For separations 1 < d < 4h^^Mpc, the change of 
peculiar velocity with redshift is only 10% of the change 
in Hubble flow. This implies a weak dependence of V12 
on cosmological expansion on those scales. 



C. Cosmological implications 



tion, as a cosmological regime, since we might be able to 
use these pairs as cosmological tracers. Measuring galaxy 
pairs in the cosmological regime would still induce a sys- 
tematic error due to the fact that peculiar velocities are 
non-zero. As the correction is of the 10% level - see 
lower panel of Fig. [2] - and we might suppose to be able 
to model this at the same level, we would have a 1% sys- 
tematic correction to contend with. This claim is little 
more than a speculation at this stage; in order to fal- 
sify or confirm it, one needs to model isolated pairs in 
detail and to analyze N-body simulations with different 
underlying cosmologies. 

Pairs isolated within a radius of 4h~^Mpc are rare ob- 
jects - see Fig. [6] - and this might result in significant 
statistical error when dealing with observations. In Sec- 
tions [VBT] and \VB2 wc found that one can drastically 
increase the number of pairs while keeping the RSD cor- 
rection small by either reducing the isolation radius to 
2h^^Mpc or allowing up to 10 galaxies to be neighbors 
of the pair. 



VI. VARYING GALAXY PROPERTIES 

The main result of the previous section, illustrated in 
Fig. |2] and [Sj is that there is a regime where isolated 
galaxy pairs may be used as tracers of expansion with 
correction for RSD that is weaker than the signal to be 
measured. Such a finding relies on the ability of mea- 
suring the redshift of galaxies regardless of their mass 
or luminosity. This is clearly not the case when deal- 
ing with actual galaxy surveys, whose flux sensitivity 
is limited. To model such a selection bias, we need to 
investigate ways of selecting galaxies from simulations 
which mimic the selection process of surveys. In this sec- 
tion we address this issue by forming subsamples where 
galaxies are selected according to subhalo mass, stellar 
mass and magnitude. We also study the redshift depen- 
dence of our results by analyzing 4 different redshifts: 
z = 0, 0.5085, 0.989, 1.504. 



Fig. [2] shows that for 0.4 < d < 4h^^Mpc isolated 
galaxy pairs at z = are nearly comoving with the Hub- 
ble flow. Thus, they move with the cosmological ex- 
pansion and, for this cosmology and epoch, only need 
a small RSD correction. Using different redshift slices as 
a way to test different cosmological expansion rates, in 
the lower panel of Fig. [5j we identify a second regime for 
1 < d < 4h~^Mpc where the peculiar velocity Vi2{z) de- 
pends only slightly on the cosmological expansion. This 
suggests the intriguing possibility of isolated pairs be- 
having in the same way on those scales regardless of the 
assumed cosmology. In particular, the lower panel of Fig. 
[5] shows that the variation in RSD model is less than the 
variation in expansion rates. Thus we conclude that there 
is cosmological signal to be extracted here. 

We refer to the intersection of these regimes, where we 
have almost comoving pairs with small redshift evolu- 



A. Method 

We select galaxy subsamples from the semi-analytic 
models by applying cuts on galaxy properties. We then 
study the dynamics of each subsample by applying the 
same analyses of Sections |IV| and |V] Initially we look at 



the number of dark matter particles of the subhalo 
the galaxy is irj^ We consider = 100, 500, 1000, 5000 
corresponding to masses m ~ 8.6 x 10^", 4.3 x 10^^, 8.6 x 
10", 4.3 X 10^^ h~^M0. Note that, where no cut is made, 
the number of particles in each subhalo is always Up > 



^ For reference, this is the np field of the Guo201C 
the Millennium simulation servers. 
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1.4 X 10"^ 


17.10 


17.81 


17.77 


17.49 


500 


4.3 X 10" 


-21.72 


1.59 X 10" 


3.1 X 10"^ 


3.85 


3.73 


3.45 


3.09 


1000 


8.6 X 10" 


-22.17 


2.60 X 10" 


1.6 X 10"^ 


2.01 


1.88 


1.67 


1.41 


5000 


4.3 X 10^^ 


-22.86 


5.36 X 10" 


3.4 X lO""* 


0.42 


0.35 


0.27 


0.19 



TABLE I: Cuts imposed on our galaxy sample from the Guo et al. (TU] semianalytic model. The first line corresponds 
to no cut at all; in that case, for the particle number column, we report the Guo et al. |10l resolution limit of rip = 20. 
The columns with a percentage sign denote the percentage of galaxies surviving the cut at a given redshift. The r-mag 
and TO* columns refer to the cuts performed at z = 0.989, as it is the only redshift we plot for these quantities (see 



Figures 11 and 12). Note that the r-mag cuts are intended to be upper limits, while the mass cuts are lower limits. 
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FIG. 6: Number of galaxy pairs for different mass cuts as a function of separation at redshift z — 0.989. In the right 
panel we omitted the m > 8.6 x 10^^ h~^M0 curve for the sake of readability. 



20, corresponding to 1.72 x lO^'^ h ^Mq, since this is the 
threshold that defines a bound subhalo according to Guo 
et al. HQ]. 

We choose the dark matter mass as our main cut be- 
cause it is directly related to the pairwise velocity dy- 
namics, which is the main subject of this paper. In order 
to make a more direct link with observations, we also 
study the pairwise velocity statistics when varying the 
rest-frame r-band magnitude]^ and stellar mass of galax- 
ies. 

The limits on r-band magnitude (hereafter r-mag) and 
stellar mass (hereafter to^,) are chosen such that, for a 
given Up cut, the corresponding r-mag and to^, cuts yield 
the same number of surviving galaxies. Table [T] reports 
the values of the limits used, together with the result- 



^ More precisely, this is the rest-frame total absolute magnitude 
in the SDSS r-band, corresponding to the r_mag field in the 
Guo2010a database of the Garching mirror and to the r^SDSS 
field in the Font2008a database of the Durham mirror. 



ing fraction of surviving galaxies at each redshift. We 
apply these cuts to the data sets before running the pair- 
finder algorithms. Thus, a pair that is isolated within 
its subsample may not be isolated when considering the 
full catalog, i. e. our isolation criterion is sample depen- 
dent. This implementation is in line with an analysis of 
an actual galaxy survey, limited by these cuts. 

In Fig. [6] we show how many galaxy pairs we find after 
imposing the cuts given in Table|l] The number of uncon- 
strained pairs (left panel) decreases monotonically with 
increasing mass cut. Note that the drop-off in the num- 
ber of pairs at large separations is due to the size of the 
individual boxes we consider and has no physical mean- 
ing. For isolated pairs (right panel), the number density 
increases as we increase the mass cut, as a higher mass 
cut results in a sparser distribution of galaxies where it is 
easier to find isolated pairs. Only for our most stringent 
mass cut, that is for to > 4.3 x 10^^ h~"'^M0, do we see a 
slight decrease in the density of isolated pairs due to the 
small number of high mass galaxies. 
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FIG. 7: Variation of the average pairwise velocity of all galaxy pairs with redshift for different mass cuts (number of 

dark matter particles) as a function of separation. 



B. Results 

1. All galaxy pairs 

In Fig. [T] we present the average pairwise velocity 
Vi2 for non-isolated pairs above different subhalo mass 
thresholds. The four panels show the same selection pro- 
cedure at different redshifts. The mass cuts range from 
1.72 X lOi°h-iM0 up to 4.3 x IOI^j^-im^ ^s tabulated 
in Table lU Note that the lowest mass cut corresponds to 
the smallest subhalo in the Guo et al. [10] semianalytic 
model, consisting of 20 dark matter particles. 

The imposition of a mass cut has a significant im- 



pact on non-linear scales. Independent of redshift, Fig. 
[T] shows that massive galaxy pairs experience an in- 
fall regime for separations smaller than d < 3h~^Mpc. 
While in such a regime, peculiar velocity increases with 
mass, with the most massive galaxies ranging from vi2 — 
330km/s at 2 = and Vi2 — 600km/s at z = 1.504. 
Lower mass galaxies, on the other hand, seem to follow 
the static solution up to higher separations, especially at 
low redshift. Our interpretation for such behavior is that 
galaxies in high mass pairs are more affected at small sep- 
arations by their mutual attraction than the underlying 
density field. 

For separations d > lOh^^Mpc, the velocity curves 
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FIG. 8: Variation of the average pairwise velocity of isolated galaxy pairs with redshift for different mass cuts (number 
of dark matter particles) as a function of separation. The isolation radius is taken to be 4h"^Mpc for each member 

of the pair. 



at each redshift seem to converge to a common asymp- 
tote. In Sec. |IV[ we have shown that this limit is cor- 
rectly predicted by linear theory - see the agreement be- 
tween the green curve and Vi2 in Fig. [l] This means 
that, even though the non- linear dynamics of the differ- 
ent mass limit pairs differs, their behavior at large sepa- 
rations seems to be predicted by the same linear theory. 

A remarkable feature of Fig. [7] is the different redshift 
dependence of the various Vi2 curves. The pairwise ve- 
locity of the heaviest galaxies (cyan triangles) greatly in- 
creases with redshift, while in the uncut case (red circles) 
it decreases. The intermediate curves seem to experience 



smaller variations. 



2. Isolated pairs 

Having analyzed the dynamics of non-isolated pairs 
with varying subhalo mass, we now do the same for pairs 
isolated within a 4h~^Mpc radius. In Fig. [s] we show 
the average pairwise velocity Vi2 for different mass cuts, 
with the redshift varying form panel to panel. This is the 
same setup as in Fig. [7] note, however, that here we only 
plot separations up to 4h~^Mpc. 
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FIG. 9: Same plot as in Fig. 8c where both the separa- 
tion and the average pairwise velocity for each curve have 
been scaled by a factor (Afrcf/M)^/'^ where M^^f is the 
minimum mass of a subhalo: Mref = 1.72 x 10^'^ Mq/H. 
Note that the errors for each curve have been scaled ac- 
cordingly. 



All curves in Fig. [8) regardless of redshift, present the 
same features found in the uncut sample shown in Fig. 
[2] and explained in Sec. |V] Namely, we see a virialized 
region on the smallest scales, an infalling regime on inter- 
mediate scales, and a roughly logarithmic growth due to 
the void effect on the largest separations analyzed. The 
main difference from the uncut case is the separations 
at which these different regimes hold. Most importantly, 
we can see that the logarithmic growth of Vi2 begins at 
larger separations for higher masses. This is intuitive, 
since we expect the mutual attraction to be stronger in 
heavier pairs, thus overcoming the void effect even when 
the galaxies are closer to the edge of the void. 

As a result of this stronger mutual attraction, the pecu- 
liar velocity contribution increases as we consider heavier 
pairs. This means that when it comes to isolated galaxy 
pairs, low-mass pairs trace the cosmological expansion 
better than high-mass ones. More quantitatively, the 
scale do where the peculiar velocity vanishes is reached 
at larger separations for massive pairs. At 2; = 0.989, 
do ranges from O.Sh^^Mpc in the uncut case to almost 
4h"iMpc for pairs with to > 4.3 x lO"h"iM0. For 
higher masses, V12 does not even cross the zero line. 

To illustrate the redshift dependence of the peculiar 
velocity in more detail, in Fig.[lO]we plot V12 for a given 
mass-cut at four different redshifts, with the mass-cuts 
varying across the panels. Increasing the mass-cut makes 
the redshift evolution of V12 more evident. As a result, 
for m > 4.3 x 10^"'^ h~^M0, we cannot identify a cosmo- 
logical regime where the pairs are comoving and have a 
redshift independent peculiar velocity. Where "indepen- 



dence" here means that the evolution is significantly less 
than the change in expansion rate. 

A closer look at Fig. [8] shows a pattern in the differ- 
ent mass cuts. We notice that, as we increase the mass 
cut, the absolute value of the pairwise velocity increases 
and the minimum shifts to the right. This suggests that 
applying a mass dependent scaling to both separation 
and velocity may stack the curves. This is shown in Fig. 
[9j where we have have taken as an example the plot at 
redshift 2 = 0.989 (Fig. 8c) and scaled both the x and 
y axis by a factor (Mi-cf/M)^/"^ where M^ef is the mini- 
mum mass of a subhalo: M^cf = 1.72 x IO^'^Mq/H. The 
physical motivation for this scaling is to have the same 
orbital period for all curves in the Keplerian regime i. e. 
on small scales. Indeed we see on this plot that such 
scaling collapses the curves specially in the infalling and 
virialized regions. 



3. Magnitude & stellar mass 

We now make a more direct link with observations and 
study the dependence of V12 on r-band absolute magni- 
tude and stellar mass. 

In the left panel of Fig. [TT] we show the average pair- 
wise velocity V12 of all galaxy pairs at z = 0.989 for 
the r-mag cuts given in Table [ij Although the curves 
retain their qualitative shape, there are two major dif- 
ferences with respect to the mass-cut sample in Fig. [7c] 
Firstly, the r-mag selected pairs have smaller average ve- 
locities. Secondly, the velocity minima are all approx- 
imately aligned at the same scale of rmin ~ 2h~^Mpc, 
while for the subhalo mass cuts the different velocity 
curves have their minima at different separations. These 
two differences are also seen when we apply the cuts in 
stellar mass (Fig. |12a ). 

The velocity differences can be explained by the fact 
that, although the subsamples chosen based on limits in 
r-band magnitude and stellar mass preserve the num- 
ber density of galaxies selected, these are not the same 
galaxies as the ones selected by the subhalo mass cuts. 
In particular, most massive galaxies do not necessarily 
coincide with the most luminous ones. In general, dark 
matter haloes trace the velocity of galaxies more directly 
than stellar mass or r-band magnitude. Cuts based on 
stellar mass or luminosity add an additional dispersion, 
affecting the position of the minima with respect to sub- 
halo mass cuts. 

The right panels of Fig. [TT] and Fig. [12] show the av- 
erage pairwise velocity V12 for pairs with an isolation ra- 
dius of 4h~^Mpc for cuts in r-band magnitude and stel- 
lar mass respectively. These plots should be compared 
with the corresponding cuts in subhalo mass at redshift 
z = 0.989 (Fig. 8c). Even though we again appreciate 
that the pairwise velocities in the r-mag and stellar mass 
cut plots are smaller, the general dynamics shown on the 
plots are the same. It is worth noting that the almost 



comoving regime mentioned in Sec. VC for each curve 
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FIG. 10: Variation of the average pairwise velocity of isolated galaxy pairs with redshift for different mass cuts 
(number of dark matter particles) as a function of separation. The isolation radius is taken to be 4h~^Mpc for each 
member of the pair. Note that panel Fig. 10a is equivalent to the upper panel of Fig. Isj 



remains unchanged both for the r-band magnitude and 
the stellar mass cuts. It is clear that the effects of galaxy 
selection (be it subhalo mass, stellar mass or r-band mag- 
nitude) play an important role in the behavior of pairwise 
velocities for isolated galaxy pairs. 



VII. COMPARISON OF TWO CATALOGS 

The results presented in the previous sections were 
based on the semi-analytic model of Guo et al. [TU]. To 
check the robustness of these results, we also compute 



the average peculiar velocities V12 for the semi-analytic 
model in Font et al. [11 . This catalog is an improvement 
over the one presented in Bower et al. [17] to better match 
the colors of satellite galaxies observed in the SDSS sam- 
ple. In order to do this, the main modification introduced 
in Font et al. [TT] is the stripping of hot gaseous haloes 
of satellite galaxies into the GALFORM semi-analytic 
model for galaxy formation, while Guo et al. [10] concen- 
trate more on the independence of satellite galaxies from 
the FOF group. Both semi-analytic models have similar 
galaxy luminosity functions that fit the data well. 

In Fig. [13] we show the average pairwise velocity V12 
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FIG. 11: Average pairwise velocity of galaxy pairs at z = 0.989 varying r-band magnitude. 




as a function of separation for all the galaxy pairs (left 
panel) and for isolated pairs with an isolation radius of 
4h~^Mpc (right panel). Each curve corresponds to one 
of the r-mag cuts in Table |l] and should be compared 
with the matching curve for the Guo et al. [TU] catalog 
in Fig. [Tl] Note that we omitted to plot vi2 for our 
most stringent cut of r-mag < —22.86 because of poor 
statistics. In general, we found that Font et al. [11] has 
significantly less bright galaxies with r-mag < —22.17 
than Guo et al. [TU], as can be seen by the large error 
bars in Fig. [13 



A direct comparison between Figures [13] and [TT] shows 
that galaxy pairs have very similar dynamics regardless 
of the semi-analytic model used. Not only do we see al- 
most the same vi2 range, but also the almost comoving 
regime introduced in Sec. jV Cj is found approximately in 
the same range. Such findings suggest that the parame- 
ters of the semi-analytic model used do not significantly 
affect the average pairwise velocities of galaxies. 
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FIG. 13: Average pairwise velocity of galaxy pairs for the Font et al. [TT] semi-analytic model aX z = 0.989. This 
figure should be compared with Fig. 11 where we plotted the same curves obtained for the Guo et al. fnT model. 



VIII. CONCLUSIONS 

We have investigated the pairwise velocities of galaxies 
with a view towards using them as cosmological tracers 
by means of an AP style test, as recently proposed by 
Marinoni and Buzzi [8] . We have analyzed the dynamics 
of such objects within the semi-analytic models of Guo 
et al. [10^ and Font et al. ^llj applied to the Millennium 
simulation [T], and studied the dependence of their rela- 
tive velocity on local density, redshift, mass of the hosting 
subhaloes, r-band magnitude and stellar mass. 

We have first analyzed the dynamics of all galaxy pairs 
at redshift z — Q (see Fig. [T]). We have found that, on 
scales d > lOh^^Mpc, the peculiar velocity is correctly 
predicted by linear theory (TSl [H]. On the other hand, 
for separations d < 3h~^Mpc, the pairs are decoupled 
from the Hubble flow and close to the static solution. 
We argue that pairs in this regime cannot be used as 
cosmological tracers (see Appendix [A]) . 

Being interested in investigating the claims by Mari- 
noni and Buzzi [H], we have studied the dynamics 
of galaxy pairs that are isolated within a radius of 
4h^^Mpc. At z = 0, isolated galaxy pairs are almost 
comoving already for separations of 0.4 < d < 4h~^Mpc 
and only need up to 20% RSD correction (see Fig. [2]). By 
analyzing redshift slices up to z = 1.504, we have found 
that the peculiar velocities are only weakly dependent on 
the cosmological expansion (< 10% variation) for separa- 
tions of 1 < d < 4h^^Mpc (see Fig. [5|. Since expansion 
is the main property characterizing a cosmological model, 
we might assume that in this regime the dynamics of iso- 
lated pairs are independent of the underlying cosmology. 
Hence, we argue that isolated pairs in this regime could 



possibly be used as cosmological tracers with minimal 
RSD corrections. 

Imposing an isolation criterion of 4h^^Mpc, as done 
in Ref. |5, greatly reduces the number of pairs (see Fig. 
|6|. We have found that one can drastically increase 
the statistics while keeping the RSD corrections small 
by either reducing the isolation radius to 2h~-'^Mpc or 
allowing up to 10 galaxies to be neighbors of the pair. 
When dealing with observations, these adjustments may 
be helpful to reduce possibly large statistical uncertain- 
ties. 

As galaxy surveys are flux limited, we have studied 
the feasibility of a measurement by varying the follow- 
ing properties of galaxy pairs: mass of the subhaloes 
that host the galaxies, r-band absolute magnitude in 
the rest-frame, and stellar mass. Low-mass pairs ap- 
pear to be the best cosmological tracers, as RSD cor- 
rections increase with mass. More precisely, a nearly 
comoving regime is reached in our analysis only for sub- 
halo masses of m < 4.3 x 10^^ h^^M©, corresponding to 
r-band magnitudes of r-mag > —21.27 and stellar masses 
of m^, < 1.59 X lO^°h~^M0 (see, respectively, Figures|8j 
11 and 12 1. We have also found that the peculiar veloci- 



ties of galaxy pairs becomes more redshift-dependent as 
we increase the subhalo mass (see Fig. [5]). Therefore, we 
suggest that isolated pairs may not be adequate as cos- 
mological tracers if their mass or luminosity is above the 
given thresholds. 

Marinoni and Buzzi [S] selected isolated galaxy pairs 
from the DEEP2 galaxy survey [181 El] with comov- 
ing transverse separation r± in the range 20kpc/h - 
0.7h~^Mpc. DEEP2 galaxies are known to reside in dark 
matter haloes of approximately lO^^Mo/h [HI HO]- The 
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results of our analysis imply that such galaxy pairs do re- 
quire corrections for evolution and cosmology dependent 
RSD component, which is significant with respect to the 
evolution being measured. 

Indeed, the primary concern for observational studies 
is the extent to which RSD "corrections" need to be mod- 
eled. Cosmological measurements from Baryon Acoustic 
Oscillation and RSD measurements on large-scales are 
reaching a precision at the 2-5% level piH^ . and it 
is therefore reasonable to suppose that this is also the 
level at which we need to understand RSD corrections 
in order to make a useful contribution to the field from 
small-scale measurements. We have investigated whether 
selection based on local density can reduce the modeling 
burden, and find that low-mass, isolated galaxy pairs are 
preferred. However, even for these galaxies, the correc- 
tions depend on sample properties, and would need to 
be recalculated for each cosmological model to be tested: 
the only currently available way to do this is via numer- 
ical simulations. We conclude that observations of close- 
pairs of galaxies do show promise for AP-style cosmo- 
logical measurements, particularly for low mass, isolated 
galaxies. However, it is likely that modeling limitations 



will continue to be the limiting factor for the foreseeable 
future. 
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Appendix A: Observing bound systems 

Consider applying the AP effect for bound systems, 
such as shown by the triangle OBC in Fig. [T4j Here 
we cannot relate Az to the proper distance using the 
Hubble parameter, or if we force this, we have to consider 
a peculiar velocity that cancels the expansion. We can 
write the observed redshift width 



n dV\\ 

c 



(Al) 



If v\\ is the orbital motion only v\\ = Vorh, then Az = 
as photons from both galaxies are subject to the same 
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Collapsed region 



Collapsed region II 




FIG. 14: This figure shows schematicaUy two collapsed 
regions. Objects A and B, which are not gravitationally 
bound, will have different cosmological redshifts given by 
za,b — ^o'^TOm '^/^i while gravitationally bound ob- 
jects B and C, shown in the second collapsed region, will 
have a line of sight component of their velocity that ex- 
actly cancels the expansion so that particles within that 
region will all have the same cosmological redshift. 



cosmological expansion. 

We can calculate a variable with the units of distance 
as in [S] 



(A2) 



But using the arguments above Aa; = and, as dv\\ 
is independent of H{z), using H{z) to translate to dis- 
tance does not provide any extra cosmological informa- 
tion. Hence any information, even if from apparent ori- 
entation of pairs, is independent of H{z). 

We now consider bound systems that have broken free 
from cosmological expansion. In general, the Hubble ex- 
pansion velocity can be defined for any pair of particles in 
the Universe. With this definition, a static body has pe- 
culiar velocities that oppose and balance the expansion. 
To see this more clearly, consider an N-body simulation 
with a comoving coordinate system. In this coordinate 
system, a body that has a constant proper size would 
appear to be collapsing. In interpreting this system from 
an N-body simulation, one might consider this infalling 
velocity as a peculiar velocity, although in effect this is 
simply balancing the cosmological expansion. One there- 



fore sees that there is a general interplay between the ex- 
pansion rate and the peculiar velocities, which must be 
included in any interpretation of data. 

Interpreting this in terms of local curvature. Fig. |14| 
shows two collapsed regions being observed. In the stan- 
dard interpretation, objects B and C, which are in a col- 
lapsed system, have peculiar velocities that cancel any 
cosmological redshift between them. The infall peculiar 
velocity must therefore be 



^pec 



-H{z)rBc- 



(A3) 



A light ray sent from B to C, will experience a Doppler 
shift due to the motion of the objects towards each other 
in addition to the cosmological redshift. Assuming that 
the light ray is emitted at a wavelength Ae™, the Doppler 
shift changes this wavelength to Xdop and the observed 
wavelength at B Xobs ■ The change in the wavelength due 
to the Doppler shift, assuming velocities much smaller 
than the speed of light, is 



Xdop — 



A. 



1 



H{z)d 



(A4) 



Due to the cosmological redshifting, the light ray is then 
observed at a wavelength of 



Xobs — I 1 



H{z)d 



X, 



where we can substitute ( A4 1 and obtain 



Xobs — Ae 



(A5) 



(A6) 



This shows that if one treats the redshift difference be- 
tween two objects as including a cosmological expansion 
component, one cannot assume that the peculiar veloc- 
ity is independent of cosmology: for bound systems their 
combined effect is zero. In an alternative and equally 
valid interpretation, B and C live in a flat Minkowski 
space-time, which does not lead to a cosmological red- 
shift due to cosmological expansion. Photons only start 
to experience a cosmological redshift once they are free 
from the bound system, and subject to cosmological ex- 
pansion: photons from B and C traveling to O experience 
the same cosmological redshift. The line of sight radial 
velocity distribution is independent of H(z) due to decou- 
pling from cosmological expansion and hence, whatever 
the true expansion rate H{z), we should expect the same 
velocities from isolated bound systems, which are simply 
all behaving as if they were in Minkowski space-time. 



